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Abstract 



We demonstrate that in the region of Bjorken xsj ~ m^/(47r_F^)^ and/or xbj ~ |t|/(47rF^)^ 
the standard xPT for the pion GPDs fails and one must perform all order resummation of xPT. 
We perform such resummation in the large- A'^ limit of the 0{N + 1) extension of the chiral theory. 
Explicit resummation allows us to reveal novel phenomena - the form of the leading chiral correction 
to pion PDFs and GPDs depends on the small x asymptotic of the pion PDFs. In particular, if the 
pion PDF in the chiral limit has the Regge-like small x behaviour q{x) ~ i/x'^, the leading large 
impact parameter {b± — > oo) asymptotics of the quark distribution in the transverse plane has the 



form (m^ = 0) q{x,b^) - l/x" ln"(6i)/&2 



2(1 + LJ) 



This result is model independent and it is controlled 



completely by the all order resummed xPT developed in this paper. This asymptotic interweaves 
with small- a; behaviour of usual PDFs, hence it depends on the scale, at which the corresponding 
PDF is defined. This is a new and interesting result in which the chiral expansion meets the QCD 
evolution. 



Introduction 

The generalized parton distributions (GPDs) |T] (see Refs. [2 [3l [H |5] for recent reviews) are ditermined 
by the matrix elements between hadronic states of the well-defined QCD quark-gluon operators of the 
type : 

0(A) =g-(iAn)7+ ^(-iAn), (1) 

which is defined on the light-cone, i.e. — (the gauge link between different points is assumed). 
The dependence of GPD^j on soft momenta and/or pion mass can be controlled by Chiral Perturbation 
Theory (xPT), see Refs. [S H E El [IHl E] . However, in the case of xPT for GPDs apart from external 
soft momenta scale, one has to take into account another momentum scale in the problem - the inverse 
light-cone distance 1/A in the operator ([1]). When these two scales are of the same chiral order xPT 
requires resummation |12j . 

In Ref. [12] it was shown that the chiral expansions of the pion PDFs q{x) (isovector) and Q{x) 
(isosinglet) can be written in leading logs as: 



n>2, even 



n> 1. odd 



(5("~i)(x), 



(2) 
(3) 



where — {^1^^ j ^iiFt^ ) is the chiral expansion parameter, i^^ sa 93 MeV is the pion decay constant. The 
superscript reg denotes the regular contributions which do not contain the (5— functions, the summation 
index n corresponds to the number of loops in xPT. The presence of the derivatives of the (5— functions 



requires the reorganization of the chiral expansion because for x 
same chiral order. 



all terms in above sums are of the 



^Usual parton distributions are just forward limit of GPDs. 
^ In Ref. [12) the coefficients -Di,2,3 were computed performing 



-1, D2 



three-loop calculations with the result: 
25 



108 



(4) 
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For the case of GPDs the problem is even more obvious. It was demonstrated in Ref. [12] that GPDs 
contain singular contributions in the kinematical variable x^j / (2 ~ xbj)) ■ For example, appearance 

of such singular terms in GPDs leads to the following contributions to the leading order amplitude for 
hard exclusive processed: 

A{U)^ f dx "f^^'^^ ^ [b^Hl/b^)]\ (5) 

■^-1 ^ ^ fe=i ^ 

where we introduced small expansion parameter 6,^ = |t|/4(47rF^)2. The presence of such strong singular- 
ities ~ 1 / in region of small ^ compensates the smallness of chiral expansion parameter and therefore 
standard chiral expansion fails. Note that in the case of the nucleon target such failure of the standard 
chiral expansion is imperative, because if t '--^ O(p^) in the chiral counting than ^ ^ 0{p^) by kinematical 
constraints. We see clearly that one has to perform the resummation of singular contributions in order 
to obtain correct chiral expansion of GPDs and PDFs. 

In present paper we perform such resummation in the large- iV limit, where N is the number of 
Goldstone bosons (A^ = 3 in real world). The chiral Lagrangian for pious is equivalent to 0(4) cr-model 
because of the homeomorphism SU{2) (g) SU{2) = 0(4). The extension to arbitrary N is done by the 
generalization 0(4) 0{N + 1). The 0{N + 1) ct— model in the large N limit can be solved by the 
semiclassical methods (see Appendix B for details on 1/iV expansion). In such approach one performs 
partial summation of xPT diagrams to all orders in accordance with the large- counting rules. We 
stress that that this partial summation respects the chiral symmetry. 

Large— extension of chiral Lagrangian and resummation of xPT 
diagrams 

We start with the description of the large— iV extension of the standard xPT. The leading order chiral 
Lagrangian describing the low-energy dynamics of pions reads |13| 

C2^^tT[{d^Udf,U^)+ml{U + U^)] (6) 

where « 93 MeV. We write in the chiral Lagrangian the physical values of the pion mass and decay 
constant because their difference with the corresponding bare constants is irrelevant for the leading log 
approximation we use in the present paper. The SU (2) matrix U in Eq. ([6]) is parameterized in the form: 

U^^ia + iTT-r), a = F,^l-TTyF^ (7) 

which gives the following result for the chiral Lagrangian: 

1 ^ 

a=l 

We have obtained the Lagrangian of 0(4) cr— model. Next we extend this model replacing the group 0(4) 
with 0{N + 1), where N has a meaning of the number of Goldstone bosons: 

1 ^ 

C2= -[d^<jd^<j + df,TT''d^n'']+a F^ ml, + J2^"^a = (9) 

a=l 

Substituting the cr— field in terms of Tr-fields we obtain 

£2 = - In-d^n- + 1 + yjf^ F^ml (10) 

^ -^TT^d^n'' + N V{n^/N) (11) 



^For simplicity we consider the chiral limit m-,r = 0. 
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where we assume the sum over the repeated indexes of the pion fields and 

, 2 



Gt 



F^/N . 



(12) 
(13) 



The Lagrangian (|f 0|1 is written in the form which is convenient for the construction of the large— 
perturbation theory. Note that large— TV expansion implies that the introduced coupling G^r is of order 
N'^ with respect to N , i.e. chiral coupling Ft^ scales as Ft^ ^ \/]V in the large N limit. 

Now, we have to introduce the generalization for the operators 0^^'^\\) (|78l79p that define GPDs 
and PDFs, see Appendix A for all definitions and conventions. The isovector operator in 0(4)— case has 
the form: 



0'=(A) = -ie[abc] T{(3, a) [/3, a] * (tt" (|(a + /3)n) 18"+ vr^ - /3)n)) 
We generalize this operator to an arbitrary N as 

0[a6]^^) = - Jc-(/3,a) * AS n''{{^{a + P)\n) W+ 7r\^{a ~ f3)Xn) 
where symbol "AS" denotes the antisymmetrization with respect to 0{N + 1) indices: 



It is clear that the correspondence between the two definitions is given by 

O^iX) = ielabc] lim O^"''!. 

Let us consider now the following matrix element in the extended theory 



dX 
— f 
27r 



gab'gba' 



q{x). 



(14) 



(15) 



(16) 



(17) 



(18) 



It defines an analog of isovector PDF q{x) for the 0{N + 1) case, the isovector pion GPD is defined by 
obvious generalization of the above equation. 

The isoscalar PDF (and its generalization to GPD) for an arbitrary N is defined as: 



dX 
— f 
27r 



where the operator is defined as: 



^a(p)^ ^ ^^ab q(^) 



(19) 



0(A) = -T{fi,a) * a{xiX)i d +a[x2X) + 7r"(a;iA)i d +tt''{x2X) 



(20) 



with 



Xi 



X2 



/3 



2 ' 2 

In order to understand the relation of the 1/A^ expansion to the standard xPT we have to inspect 
the diagrams of the large— approach. To the 1/N order for GPDs one has to compute the diagrams 
shown in FiglT] The first diagram is the tree contribution to the pion PDF, the second and the third 
diagrams are the first non-trivial 1 /N correction^. The solid line denotes the usual pion propagator and 
the dashed line is the propagator of the auxiliary field ip which is suppressed by 1/N (see details of 1/iV 



*To the 1/N order there axe another tadpole-hke diagrams, however they do not contribute to the leading logs. 
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Figure 1: Diagrams which contribute to 1/N accuracy to the pion GPDs. 



expansion technique in Appendix B). The exphcit expression for the propagator of the auxihary field is 
given by 



O 



-f 



AF2 



N 



O 



A 



9 ? 



f M ' r(f + £)r^(f-£) 

£ [-p^Y r(2 - 2e) 



(2f) 
(22) 



where we apply the dimensional regularization with I? = 4 — 2£ and with renormalization scale /i. After 
expanding the above expression in powers of pion coupling N/F^: 



( ^ 



AF3 



o 



N 



2AF3 



(23) 



we observe that this expansion produces the infinite series of the bubble chains with the multiplicative 
factor proportional to N/F^. This expansion explains which diagrams of the usual xPT are taken into 
account at given order of 1/N expansion. 



Large— A/^ resummation for the pion PDFs 

In this section we consider the large- iV resummation of the singular contributions to the pion PDFs. We 
give detailed calculations for the case of the isovector PDF q{x). For the isoscalar case we outline main 
differences in the calculation comparing to the isovector case. 



Isovector PDF 

Let us compute the leading logarithms from the diagrams with the insertion of the bubbles and select 
the contributions with the derivatives of the (5— functions. Only two first diagrams in Fig. [1] contribute 
to the isovector PDF, the third diagram obviously contributes obviously to the isoscalar operator only. 
The corresponding diagrams are [/T^— divergent and must be renormalized. The t/V— subtractions is 
performed by the application of R-operatiorif|. The series of the relevant diagrams originating from the 
loop diagram in Fig[T]can be written as 



\2 2 



O 



(24) 



The (5— function in the numerator arises from the non-local vertex of the operator. It is clear that each 
term in the sum (|24p corresponds to n -I- 1-loop diagram. The R-operation removes the C/y-divergencies 
in the n + 1-loop diagram. The subtractions have following structure 

n 

R g("+i) = G("+i) + \G^r.+i-k) ^ J2 ^cL'i , (25) 

{Ffc} ^ k=l ^ 

where G^"^^-* denotes the bare (n -t- l)-loop diagram, the second term describes the subtractions of the 
subdivergencies in the bare diagram, the sum in the second term runs over the set of the C/y— divergent 



^An introduction to tlic renormalization in quantum field tlicory and detailed discussion of the R-operation see in |14II15| . 
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Figure 2: The diagram with the n— loop t/]/— counterterm. Corresponding counterterm is shouwn by 
the dashed Hne 



fc-loop subdiagrams {Ffc}. The index k denotes the number of the loops in the subdiagram and the 
pole originates from its C/V^— divergence. It is clear that the index k takes values 1, . . . ,n. The term 
Q{n+i-k) (^QgQYibes the (n+1 — fc)-loops diagram with insertion of the counterterm vertex that results from 
the contraction of the fc-loop subdiagrams. The last term denotes the subtraction of the total resulting 

(k) 

divergency, G^j-gg is the notation for the tree level diagram with the appropriate vertex structure. The 
contributions to the coefficient in front of the leading logarithm arises from the terms in (|25p except the 
last one. The e pole structure of the (n + 1 — fc)-loop diagram G''-"+^~'^^ in Eq. is the following: 

= ^ (^) G,ie) (26) 



^ (n + 1)! 

In the above equation we took into account that Gk{s) is the polynomial of order (n + 1 — fc) in e and 
it is an analyitical function of external momenta and the pion mass. In the sum of all contributions in 
Eq. (|25p the poles in Eqs. with non-analytic dependence on external momenta and rriTr must cancel 

(k) 

due to locality of the total counterterm. That imposes certain relations between the coefficients 
with k = 0, . . . , n and / = 0, . . . , n. One can easily show that these relations leave only one unknown 
coefficient zhat must be computed in order to find the total coefficient in front of ln'^"^^'(/i2/m^), and 
has the form: 



9 ; / 9a - — 7 TT, + • ■ ■ ■ (27) 

It is convenient to introduce coefficient gQ as an unknown quantity, i.e. one has to compute only the 
contributions X){r } p^^^^^- -'^^ other words, we have to compute only one loop diagram, but with all 
possible n— loop counterterms. An example of one of such diagrams is shown in Fig. [2l The structure of 
the n— bubble counterterm has the form: 

^ {lnikpr+' + + 7(°'^)P] (M" + • ■ • ) (28) 

where the dots denote all other possible contributions constructed from the m^, k^ and {kp). The power 
n is dictated by dimensions. Inserting the expression (j28p into the original diagram instead of the bubbles 
chain one obtains 

G(^) ^ F'-\(3, a)^l£^ 'lk2-%~/''^^ {inikpr^' + [iL'^V^ + [kpT + . . . ) (29) 

Now we observe that the most singular term (with the highest derivative of (5-function) is produced from 
the contribution with the highest power of the scalar product {kp). Expanding the (5— function in the 
numerator of Eq. (p9| 

k+5{l3k+ - xp+) = kl+^ f3''6^''\xp+) + ... (30) 
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Figure 3: The diagram with the n— loop f/y— counterterm which includes the operator vertex and 
resulting reduced one-loop diagram. The counterterm is shown by the dashed hue 



and neglecting the terms with lower power of fc+ , one obtains 

Computing the one-loop integral over k we obtain the expression for ^q"-* = (n + 1)! 7„. Using Eq. (P?]) 
we obtain the desired singular contribution to PDF: 

q{x) ^ 7n<5(") {x)a'^+' ln"+i [1/aJ + • ■ ■ , (32) 

where we fixed the renormalization scale /i — [AttFt^)'^ for simplicity. 

Going back to the general case, we compute the sum over all possible n— loop subdiagrams. For that, 
in addition to the considered case, one has to consider the subdiagrams with the operator vertex, for 
instance such as in Fig. [31 Obviously, the reduced diagrams have always at least one contracted line, 
this reduces possible power of derivatives of the i5— functions in the resulting expression. Indeed, the 
contracted propagator is always proportional at least to one power of the pion mass that excludes 
one power in possible expansion of the (5— function. In other words, there is no possibility to produce the 
singular contributions of type (I32p from the contributions with vertex subdiagrams in the diagrams with 
the bubble chains. 

From the above consideration one learns several important lessons. We expect that each term in the 
sum of diagrams p4)) contains a singular contribution proportional to (J*-"-* {x)a^'^^ In""*"^ [1 / a^] . Moreover, 
taking into account the relations due to the cancellation of the poles with non-local dependence on external 
momenta, we obtained that the calculations of these contributions can be reduced to a consideration of 
the diagram as in Fig[2]with rt— loop bubble chain counterterm. Important technical observation, is that 
the structure of the corresponding counterterm is restricted by the maximal power of the scalar product 
(kp) as in (j29p and does not depend on the pion mass in to the subdiagram. Therefore in the practical 
calculations we can put the pion mass inside the bubble chain to zero. In addition, one can expand the 
(5— function in expression ()24|) and pick up only the contribution with (5'^"^(x). Using these simplification 
and denoting the most singular term as q™'^{x) one obtains 



(n+l) / ^ n+1 



5sing(^) ^ 1 y (X")^(")(xp+)R / — ± ^ ( ^[{kp) 



o 



(33) 

J m^=0 



where the R-operation removes subdivergencies only from the bubble chain subdiagram and 

pi 

(x") = F'=\(3, a) *(3'' ^ J dxx" q (x) . (34) 

It turns out that such simplified expressions can be computed directly, without the combinatorial analysis 
of the R-operation mentioned above. We provide the corresponding technical details in the Appendix B. 
The result is very simple 

r, °° n+l 

n=0,2,4,... ^ 

where e = ^-a^ Inll/a^]. This expression shows that one has the singular contributions to all orders and 
we can perform the summation of these terms. Let us represent the n-th derivative of the i5-function as: 

= f ^(jA)"e'^^ (36) 
7 2n 
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Changing the order of the integration the summation we obtain 



-^-^ n=0,2,4,... ^ ' 



E =r^-E^^=-r^--^- (38) 



-1 



Computing the sum as 

oo 

^ (n+ 1)! y 1 ^ n\ 2 

and inserting it into (j37p we obtain 

11 1 ^ 

^.ing(^)^_i^ /• ^^"(^^Z' dTd{x + PeT) = -^ f dp^9{\x\/e<m. (39) 



Taking into account that isovector PDF is symmetrical function q{—x) = q{x) one can rewrite ([39]) as 



This is our final result for large- TV resummation of chiral corrections to the isovector pion PDF. The above 
equation demonstrates explicitly that resummation solves the problem of the singular terms - the infinite 
sum of (5-function derivatives results in the contribution of the type f{x/a^)/a^ as it was conjectured in 
Ref. [12]. Note that the large- resummation preserves all chiral properties of the considered quantities. 

From l[40|l . one can observe that resummation does not change the small- a; asymptotic of the PDF. 
Taking the Regge-like small- a; behaviour of PDF in the chiral limit: 

^ (-) - ^ (41) 

we obtain non-analytic behavior of PDFs in the chiral coupling (a^ ln(l/ax))'^, for the chiral correction 
to PDF in the small— a; region: 

,sing(^) ^ r ^ "q^^^^f «xln(l/«x) r (42) 



Interestingly, the leading power of chiral coupling is determined by the intercept of the Regge trajectory 
uj {u! Ki 1/2 for the isovectors case). This shows clearly importance of the resummation of singular chiral 
corrections for the derivation of the leading chiral counting of PDF. Naive chiral counting - without taking 
into account the second scale related to the light-cone distance-suggests that the leading chiral correction 
to the isovector PDF is ^ a^ln{l/ax)- However, it is not correct, as the leading chiral corrections are 
dominated by the small- a; region and one must perform the resummation discussed in the present paper. 

Isoscalar PDF 

In the case of isoscalar (C parity even) pion PDF one has to compute additionally the third diagram 
shown in Fig. [TJ The dashed line in this diagram represents the propagator of the auxiliary field (see 
Eq. I2ip corresponding to the sum of bubble chains. The third diagram is more complicated as it contains 
two loop integral. The expression for this diagram has the form: 



N ^ 

n,m— 

1 



TVs m+n+2 



k+6{k+P - xp+)[P - m2]"+'"+2 



r(i + £)r2(i-£) 
r(2-2£) 



m+n 



(43) 



(44) 
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P-K/2 > i' < (X) *' • > P+A/2 

k+A/2 ((-A/2 



Figure 4: The diagram for correction for isovector GPD with momentum flow. 



Applying the renormalization procedure (R operation) and picking up the most singular in x terms we 
obtain the contribution of the third diagram to the singular part of the isoscalar PDF Q^"^^{x): 



_o °° i^a ln(-!-))"+"+^ 

- E '-"'">"'"'""W (,„:„^3)T(„ + ,„ + i) • <•" + " + '> - <«> 

n,m— 

Each term in the above sum depends only on m + n combination of the summation indice; hence, this sum 
can be easily computed. The final result of the summation of both the second and the third diagrams in 
Fig. [T] is the following: 

Here, as before, e — -j'^x ^'^(l/'^x)- ^^'^ '^^ isovector PDF, the small x behaviour of the chiral 

correction has the same form as for the PDF in the chiral limit, e.g. the Regge -like that ~ 1/x'^. One 
can easily calculate that the small-x asymptotic of the leading chiral correction has the form: 



y ^^ j ■ (47) 

Again we see that the correct dependence of the small- a: PDF on the pion mass can be obtained only 
after the resummation of all orders of the standard xPT. 

Eq. (|42p and Eq. (|47p show that the leading chiral dependence of the PDFs on is determined by 
the small x behaviour of the PDFs. The dynamics of the chiral degrees of freedom (Goldstone bosons) 
and the dynamics of QCD degrees of freedom (quarks and gluons) intertwine. Noticeably, even the power 
of the chiral expansion parameter changes with QCD evolution of PDFs. 



Large— resummation for the pion GPD 

In this section we consider the resummation as before but for the GPD matrix element. The diagrams 
describing the 1/N correction to the pion GPDs H{x,^,fB are the same as in Fig. [H bu with momentum 
flow depicted in FigHl The expression for the 1/iV correction to the isovector GPD has the following 
form: 

oo 

H'=\x,tt) = -Y.KF'-\f],a)* (48) 

f d^k k+S{j3k+ + < - xp+) f N , y+' r 1" 

J {27:)D [(fc- A/2)2-m2][(^ + A/2)2-m2] {f^^^ ^ ^' ^ ) V \ ' 

Due to the presence of the momentum transfer instead of (5-functions we obtain the contributions ~ 
^^2J<iI, Such contributions in the forward limit (i^ — > 0) are reduced to the derivatives of (5— function. 
Therefore we shall separate them neglecting the others less singular terms. The explicit resummation 
can be done in the same way as for the PDF. But now the calculations are slightly more complicated 



^We do not introduce special notation for this correction in order to avoid proliferation of notations. It is clear from the 
context which part of H{x,i;,t) is considered 
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due to additional parameters ^ and t. At the end we obtain the following expression for the singular 
contributions to the isovector pion GPD: 



where 



ml - t{l - 7f) / 4, 



R = 



R\t=o = R\jn^=o = (1 - "rf ) b 



X' 



is the small chiral expansion parameter. One can easily perform summation of the singular contributions 
in Eq. (|49|) with the following result: 



H'=\x,^,t)^-^ j'^dri j^^dr (^l + !^jF^=i(;3,a)*J([?77 + r]/3 + ^a-x), (50) 
where 

C = yi?ln(l/i?). 

The result of the large -N resummation of the singular contributions to the isoscalar pion GPD has the 
form: 

il^=°(x,e,i) = J ^dv dr {3-2T/0 F'=°{P,a)^f^(^l+^y{[7j^ + r]P + ^a-x){51) 

I-—] S {[t]^ ~T]/3 + ^a-x)~ (77^/3 + ~ x) 

T ) T 

Note that the small chiral expansion parameter e = ^Rh\(\/R) enters these expressions always divided 
by ^, meaning that for ^ ~ e the singular contributions are not suppressed by the small chiral parameter. 
In next sections we shall analyze obtained expression for resummed chiral corrections to GPDs and 
corresponding amplitude of hard exclusive processes on the pion target. 

Parton distribution in the transverse plane 

To analyze the distribution of partons in the transverse plane one has to take the limit ^ — > keeping 
< 7^ in GPDs. One can easily perform this limit for the expressions (I50l5ip with the result^: 

^ dp o 



H'=\x,^ = Q,t = -Al)^q^^^{x,A^)=-^e{\x\<e) ^q{P) ^1 - J, (52) 



and 



//i=O(x,e-0,t--Ai)^Q-"s(x,A^) = -±9{x<e)j'^Q{P) (53) 



'1 + Jl 



^In section we consider the case of m-n- = 0, because we focus our discuss on the dependence of GPD on t and ^. For the 
real world this case corresponds to the parametrically wide domain of momentum transfer squared <C — i <C {47rF,r)^. 
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where e = ^-h^lnil/b^) with 6^ = |t|/4(47ri^^)2 = Ai/4(47ri^^)2. 

The Fourier transform of the PDFs q{x,t), Q{x,t) in A_l (A^ = — t) gives the density of partons in 
transverse plane {b±) with given momentum fraction x [16j . The chiral expansion allows us to compute 
the asymptotics of the distribution in the transverse plane at 6_l ^ oo. The one-loop xPT gives the 
following large b± behaviour of isovector quark density q{x, b±^j^ [7]: 

qix,b^) ^ ^ Six) . (54) 

The large b± behaviour of isoscalar distribution of quarks Q{x, b±) in two-loop xPT has the form: 



g(x,6^)^-^5'(x). (55) 

We see that at large b± the density of (anti)quarks is formed by partons with small longitudinal momentum 
fraction x. The presence of the ^-functions is the artifact of the finite order xPT calculations, we simply 
can not resolve the 5- functions. This shows that we must perform the resummation of xPT in order to 
access the correct large b± asymptotic of quark density in the transverse plane. We can easily perform 
the Fourier transormation of our results for q{x,t) and Q{x,t), see Eqs. (|52l53p . As we discussed above 
only partons with small x are relevant for the large b± asymptotic of parton densities. Assuming the 
Regge-like behaviour for the forward parton density q{x) ~ and Q{x) ~ (the subscript ± 

stays for the C-parity [signature] of the Regge trajectory) we obtain the following result for b± oo: 



^(-^M-||yS^^, (56) 

' ^ {b\y+'^+ x'^+ 

The intercepts of the Regge trajectories are uj- 1/2 and a;+ w 1.1. We see that distribution of partons 
at large b± (|56p drops slower than the naive results (|54l55p of finite order xPT. Also we observe very 
interesting phenomenon - the distributions of partons in the transverse plane at large impact parameter 
b± depend on the intercept of the corresponding Regge trajectory. This new phenomenon is revealed due 
to the all order resummation of xPT developed in this paper. 

The average transverse square size of a hadron in the transverse plane is defined as: 



(bDix) = J d^bi_ bl qix,b^). (57) 

For the isoscalar PDF Q{x,b±) the integral for {b']_){x) is IR- finite in the chiral limit (see Eg. [55|l . meaning 
that the corresponding radius is not determined by the chiral expansion. On contrary the integral for 
the isovector {b\){x) is IR-divergent and, hence, its value is dominated by the large b± ^ 1/m^, so that 
we can obtain the behaviour of {b']_){x) in the isovector case: 



We see again that at small x the transverse size of the pion grows as the powei0of 1/x, also its dependence 
on small pion mass is power like (compare this with ~ ln(l/m^) behaviour of the isovector radius of the 
pion in the standard xPT). The growth of t he nucleon transverse size at small x has been discussed 



*We remind that we consider the chiral hmit m.^ = 0. 

^Note that in the Regge-like Ansatz for the t-dependence of GPDs the transverse size of a hadron grows as ~ ln(l/a;) at 
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in Ref. [17]. The power- like growth of the pion transverse size at small x derived here implies that the 
similar phenomenon takes place for the nucleon case0. 

Up to now we discussed the large 6_l asymptotic for the chiral limit TOtt = 0. The case of 7^ can 
be obtained from from the general result given in Eqs. (|50|5ip : 



TV In [bljAixF^f] 



(59) 



For rriT^bj^ 



( h ^ 1 1 [ N \n[bl{4nF^)^] 
1^'^'''^)- Nbl[\x\ bliAnF^r 



0{m^b^). 



This result coincides with the asymptotics (|56p for the massless pion. For m„6 



q{x,bj_) 



1 e 



N 



N m^b± In [b]_{4:TTF^y 
R 5i(47rF,)2 



1 + 



m-^b^ 



Chiral corrections to the amphtudes of hard exclusive processes on the pion 

In this section we compute the contribution of the singular terms to the amplitude of hard exclusive 
processes. We consider the Mandelstam t in the region ^ —t ^ (47rF^)^ in order to focus on the the 
i-dependence of the amplitude, i.e. we can put effectively ~ 0. 

The corresponding amplitude is computed in terms of GPDs as the following convolution integral: 



A'=''\U) = ^^^^-^^ ■ (60) 
J^i £.~x-iO 

Let us start with the isovector amplitude. To calculate the contribution of the singular terms ~ 
[6;^ ln(l/6^)]'^/^A: into the amplitude we substitute the expansion given by Eq. ([49]) into the integral 
([60]) . The result can be represented as the following series: 



N ^ \ ^ \by (n+1) £"+1 7 1 \-z 

n=0.2 ^ \ X/ / \ < J S J-L 



(61) 



Here we introduced the distribution amplitude (DA) of two pions with angular momentum I = n + I 
defined as: 

where Pn+iirj) are Legendre polynomials and the two pion DA ^18J $ (2,77) is expressed in terms of 
coefficient function F{/3, a) as follows: 



oI=l 

$ {z,ri)^T]F^=\a,P)*S{/3ri + a-z). (63) 

The sum (j6ip has the form of the partial wave expansion of the amplitude in the t-channel, see |19) . 
Each term of the sum (j6ip is real, however for £ < Nb^ In (1/6^^) the sum is divergent and the imaginary 
part of the amplitude is generated - the case typical for the duality. That naturally invites one to apply 
the methods of the resummation of the ^-channel exchanges developed for the dual parametrization of 

^"^We note that the results obtained here for the isoscalar quark distribution are also valid for the gluon distributions. 
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GPD in Ref. [20] . The final result for the chiral expansion of the imaginary part of the amplitude is the 
following: 

W-^(^,t) . . .(^ < 2.) J^'^^^ '-^ N (X) (l - , (64) 

Here e = ^b^\n{l/b^) [b^ = |t|/4(47r^V)2], TV (a::) is the so-called GPD quintessence function [2H [22]. 
which is expressed via the imaginary part of the amplitude at = i = with help of Abel tomography 
methods Ell [53]: 

Im 1(0 = TT if (0^. (67) 

o 

Given the Regge-like small ^ behaviour of the amplitude at t = -Im A {0 ^ ^/^'^- The leading chiral 
correction to the amplitude has the form [12] : 



Im AiC,t)^ — 



(68) 



which depends on the intercept of the corresponding Regge trajectory (w w 1/2 for I ~ 1 and a; w 1.1 for 
1 = 0). This novel phenomenon can be obtained only after all order resummation of xPT. The standard 

xPT gives the leading chiral correction of the form Im A^^^{^, t) ^ b^ In (^y'^ and Im ^^^°(0 t) ^ 



K In 



Conclusions and discussion 

We demonstrated that in the region of Bjorken xbj ~ / {AttF-,^)'^ and/or xbj ^ \t\/ [AttF-j^Y the ordinary 
chiral expansion for the pion PDFs and GPDs fails and one must perform all order resummation of 
the standard xPT. We perform such resummation in the large-A^ limit. Explicit resummation allowed 
us to reveal novel phenomena in quark mass expansion of PDFs, in low energy behaviour of GPDs 
and amplitudes of hard exclusive processes. The main qualitative (model independent) results are the 
following: 

• The leading small ttItt asymptotic in the region of small x of pion PDFs depends on the intercept 
(w) of the corresponding Regge traiector-\Pl: 



1 



q{x) 



(4^F,) 



■In 



1 



(69) 



• The leading large impact parameter (&j^) asymptotics of the quark distribution in the transverse 
plane has the form (we show result for m,r, the result for 7^ is given in Eq. (j59|) ): 



q{x,bi) 



l_ ln"(bi) 
(6i)i+- 



x'^ 



(70) 



^^We consider the Regge-like behaviour of PDFs for simplicity, one can easily obtain corresponding leading chiral correc- 
tions for other types of small x behiour of PDFs. 



12 



The distribution of quarks at large impact parameter is controlled completely by the all order 
resummcd xPT developed in this paper. This asymptotic is determined by the small-x behaviour 
of usual PDFs, hence this asymptotic depends on the scale, at which the corresponding PDF is 
defined. This is new and interesting result - the chiral expansion meets the QCD evolution. 

• The leading small t behaviour of the amplitude for hard exclusive processes on the pion target has 
the form: 



(71) 



Measurements of such processes at small xbj and small t would allow us to probe the chiral dynamics 
in a completely new regime - the dynamics of chiral and quark-gluon degrees of freedom intertwines. 

The complete results for the resummation of xPT for pion PDFs and GPDs are given in the main body 
of the paper. 

In order to perform the resummation of xPT we used 1/iV expansion. To assess the accuracy of the 
expansion we computed the first three coefficients in Eq. ^ for arbitrary TV with the result: 



For = 3 the above result coincides with the three loop calculation (see Eq. (U))) of Ref. [T^], also for 
1 the D2 coefficient is zero as it should be, because the 0(2) (T-model is a free field theory. All 
that provides strong check of the three loop calculations of Ref. |J2j. From Eq. (|7^ we see that the 
corrections are sizable (30-40 %) and seem to increase with number of loops. To check this tendency we 
computed four and six loop D4 and Dq coefficients in the expansion ^ with the results:) 



76 1 113 1 64 1 \ 3 
27 '2nlP " 27 J 

46208 1 9441001 1 36581227 1 30246239 1 2449121 1 \ , 5v 
10125 iV ^ 911250 7V2 ~ 2733750 ^ 2733750 ~ 546750 J ' ' 

These coefficients are zero at iV = 1, that provides a check of our calculations. We see that indeed 
the 1/A'^ corrections increase with number of loops (they amount to 60 — 75%). However, we note 
that the total coefficients £)„ at TV = 3 decrease rapidly with number of loops, that indicates that the 
resummation of singular contributions is possible beyond the large-iV expansion and that one expects 
that main qualitative conclusions of the present paper are valid also for the exact resummation. 



7N^ 
480 

8960 
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Appendix A. Twist-2 operators and their matrix elements 

In this Appendix we briefly describe the definitions and some technical details used in the paper. We 
introduce two light-like vectors n and n: 

Ti^ — fi^ — 0,n ■ fi — 1, a+ = a ■ n. (73) 
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There exist two QCD quark light-cone operators of twist-2: 

^fl = ^(l-75), Pl = ^(1 + 75) 



(74) 



[OR]fg = Qa (5A n) qf (-^A n) , (75) 

[OL]fg = qg (iA n) {-^X n) , (76) 

where indexes f,g stand for flavor. These operators transform under the global chiral rotations as 

Ol ^ VlOlVI , Or ^ VrOrV^ . (77) 

In xPT these QCD operators are described by an effective chiral operator with unknown chiral constants. 
In the pure pion sector one finds [3 [HI [5] 



4 



Of, (A) = -'-^T{P,a) 



U 



a + /3 
2 

a + /3 



Xnj n- d W 
Xn]n-dU 



a — (3 
2 

a — /? 



An 
An 



fg 



fa 



where by asterisk we denote the integral convolution with respect to (3 and a: 



T{(3,a)*0{f3,a) 



dp 



1-1/51 



1 + 1/31 



da T{I3, a) 0{I3, a) . 



(78) 
(79) 

(80) 



Here J^{(3, a) represents the real generating function for the tower of low-energy constants and denotes 

a combination of derivatives — d^. It is important to note that the expressions ([75]) and ([TS]) describe 
correctly only the operator vertices with two attached pions (including the tadpoles loops). The low- 
energy constants ^(/3, a) characterizes the structure of the pion, they are not determined in the effective 
field theory. According to the isospin / = 0, 1 one can construct two independent functions: 



which are convenient for description of the pion matrix elements. Pion PDFs are defined as 



277 



{'k\v) |tr [T^OL+fl(A)]| -K^v)) = ^ie{ahc\q(x) 



dX 



-ip+ a:A / b 



{n\p) |tr [OL+fl(A)]| n^{p)) = 26->'Q{x) 



which in the chiral limit can be written as 

-•1-1/31 



-1+1/31 
1-1/31 



da F'="(/3,a) 



9{P) q{p)-e{-f3) q{-(3) 



da F'-\p,a) 



-1+1/31 

The GPDs are defined as: 
dX 



0{P) Iif3)+e{-l3) q{-P) =q{f3). 



iP+ x\ I h I I 



2ti 



{n'ip') |tr [T^OL+RiX)] I 7:^{p)) - Aze[abc]H'-\x, ^, t), 



(81) 
(82) 

(83) 
(84) 

(85) 
(86) 

(87) 



dX 
— ( 
2tt 



(^^(p')|tr[Oi+H(A)]K'^(p)) 



25°''iJ^=0(a;,^,t) 



withP- i(p + p'), i 



_ (p'-p 



(p'+p)^ , t — {p' — p) . In the forward limit ^ ^ 0, t — > 0: 
H^='^{x,Q,0) = q{x), H^=°{x,0,0) = Q{x). 



(89) 
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Appendix B. Large N expansion in 0{N + 1) model 



In this Appendix we discuss the technique of expansion for the 0{N + 1) sigma model. We follow 
closely the methods of expansion presented in the book by A.N. Vasiliev [14] . 

Diagram technique for the large— expansion 

Our task is to construct the large— expansion for the field theory given by the Lagrangian 

C = -^TT''d^TT'' + NV{n^/N) + J''n'', (90) 

V{n'/N) = I ffi^) + ^"G^G^ ^ G' = F'lN^O{N') (91) 

where J° is a source for the pion field. The constants m and F represent the bare values of the pion 
mass and axial coupling. 

Introducing two auxiliary fields '0 and ip and using the representation 

1 = j D^j (5(7r2 - TV V) = j Dtjj Dip exp [i (p {N - n^) /2] 

one can rewrite the Lagrangian as 

£ = -^TT^d^TT" + p{N^j- TT^) /2 + NV{il;) + J^tt", 

Performing the integration with respect to the pion fields we obtain the effective Lagrangian 

Ceff = A^^tr In [{d^ + p>) /d^] + N ip^/2 + NV{yj) + .^{8^ + p>)-^ r /2 . 

Notice that the second operation arises from the normalization constant in the functional integral. 
Computing the resulting integral by the steepest descent method, one obtains the desired large— 
expansion. 

The equations of motion define the constant solutions ip — po and ^ = ■ 

Let us postpone the discussion of the solution and suppose simply that it exists. Then performing 
expansion around the classical values we obtain the effective action 

Ceff = nUiIh [{d^ + P0 + p)/d^] +N {p + po) + Vo) /2 (93) 
+ NV{i^ + i;o) + J''{d^ + Po + py^ J72 (94) 

= -^p{x)iA'^{x, x)p{x) + N p{x)^{x)/2 - ^ipK^^p (95) 

- r{x)A{x, y)p{y)Aiy, z) J«(z)/2 + ... (96) 

where we introduced convenient notations: 

= \{G^- ^o)"' d\ {d^ + ^o)-' = i A(x, y), 

In the second line (jOS]) we show the relevant to our accuracy 0{N^^) contributions only . The first term 
in (j95p arises from the trace in the first line, and A'^[x, x) in graphical notation corresponds to the scalar 
loop 

A\x,x) = L^ I „ ' 1- . (97) 



15 



The quadratic in ^ term originates from the expansion of the interaction NV{^) and the dots denote the 
terms that contribute to the higher order corrections ~ 0{N~'^). We also keep the contribution quadrat- 
ical with respect to the sources because we consider only the two-pion matrix elements. Performing the 
diagonalization with the help of substitution ip ^ ^ + ^K^^tp we obtain a simple expression for the 
relevant part of the effective action 



r ^ 



2 92 



ifi-J^AipA J72 + ... . (98) 



The quadratical term provides the propagator for the auxiliary field ip. Notice that insertion of this 
propagator produces the factor 1/N. This provides the counting rules with respect to N. The diagrams 
in FigUare produced when we contract all the auxiliary fields, differentiate with respect to sources and 
amputate obtained irreducible diagrams. One can easily check that the terms denoted as dots can provide 
only graphs which are suppressed as 1 /N'^ . 

Let us shortly discuss the solutions of the equations ([HH). It is clear that constant ipQ provides the 
physical mass ttItt of the pion in the large— iV expansion and the second constant tjjo can be associated with 
the physical axial coupling Ft^. On can easily see that all diagrams which contribute to the renormalization 
of the mass and coupling can not provide leading log's to the singular terms because renormalization 
decreases the power of the logarithm in front of (5— function. Therefore we can skip the discussion about 
the renormalization of the pion mass m and axial coupling F. Moreover, from the consideration in the 
text, we know that in the four-pion blob we can neglect the mass of the pion because corresponding terms 
produce subleading logarithms only. Therefore, in the propagator of the auxiliary field ip one can neglect 
the pion mass. That is equivalent to 

^o^rnl^ 0, (99) 

It is clear that we must keep the pion mass only in the propagators connecting the operator vertex with 
the four-pion subdiagram. Next, in the propagator we substitute 

(G^ - ^0) - F^/N = Gl (100) 

Finally, we obtain the following simple Lagrangian 



r ^ 



2 9^ 



ip-J'^AipA J72 , (101) 



where Lq denotes the massless scalar loop ((97)) . Thus for the propagator of the auxiliary field one has 

(<n<n) ^ — P^ 1 . (102) 

LJ^ 1 -I- -^-Ijq 

Calculation of Large— corrections to PDF 

Our aim is to discuss the evaluation of 1/A'^— corrections given by graphs in Fig[T] Consider one loop 
diagram which contributes to q^™^ in ([55]) . Expanding the propagator (|102p in series with respect to loop 
Lq we obtain expansion of q™^^{x) which can be written as 

2 °° 

g""S(x) = --^(:r")<5(")(a:) J„. (103) 
11=0 

Hence, we have to work with the expression which include insertion of the products of the n massless 
loops (^qT'^Lq^ in each J„. It is clear that such construction is associated with the n + 1— loop graph. 
The renormalization of such graph includes the substruction of various subdivergencies. To the leading 
logarithmic accuracy it is enough to consider the insertion of the Lo~counterterms only. The explicit 
expression for Lq reads 



F2 2 ° F2 2 (47r)2 ^-p2y £ mle^ ^ 



(104) 
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where for simplicity we introduced convenient notation 



2 {A^F^Y 

For the product of n loops Lq we obtain obviously 



N m? rfl + e)rfl - e)2 



P 



r(2 - 2e) 
2 ^"e^(e)"" 



The expression for the total unrenormalized graph with n loops Lq can be represented as 



^n+l 



i {-k+/p+) 



(n+l) 



(A-e) 



2{pk) 



mt 



(n+l) 



(105) 



(106) 



where A = {n + 1) e. Evaluation of this integral is standard. Rewriting the denominator with the help of 
Feynman parameters, shifting the integral momentum k and evaluating the resulting integral we obtained 



_i n" 1 (fJ-^/ml) r [1 - A, 1 + A, 15/2 + n + 1, 3n + 5 - 2A] 



e" A n\{n + 2) F [n + 2 - A, D/2, 3n + 5-A - e] 

2 X A 

En [s, A] . 



1 



(107) 
(108) 



e"+i (n+l) \ml^ 

Subtraction of Lq— subdivergencies in the n— loop chain [Lq]" is described by the substitution (|104p : 



N ip^ 
F3^ 



E 

fc=0 



9 \ n n 

p e 



E 



k=0 



TV 

F3~ 



n 



I — A; 



(109) 



,2 \ e 



-p^ 



Each term in the last sum can be represented in the form p06p with A ^ A' = {n+l — k)e. Hence after 
subtraction of the subdivergencies we obtain following expression 



fc=0 



(110) 



Substituting the J— integral as in (|108p we obtain 



Jn — 



fc=0 



1 E ( k ) ( m2 



2X^' (-1)^ 



(n - fc + 1) 



En [e, A'] . 



(Ill) 



Performing expansion of this expression with respect to e we obtain the local poles and finite terms: 



Jn — 



f + ... + ^ + ln"+i [fiym'l + 0(ln" [^Vm^J ) 



~n+ 



(112) 



In order to compute coefScient J^^ we take into account that the largest power In"''"^ arises only from 
the expansion 

,2 X A' 



1 



-n+l 



il^^A±Ipln"+i [^V-^] +0(ln" [mV-^]) 
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Picking up coefficient in front of In"^^ in we obtain: 

^n[0,0] y 
(n + 1)! ^ 



- (^In [,Vml]r' ^ ±(l) i-lt in~k + If 

fc=o ^ ^ 



(n + 1)!' 
where we used that 



1 _ k)- ( l)-n\, Ihn i?„ [e, A] = 1 (113) 



k=0 

and redefining eln [/^^/m^] = e one has 



Then the expansion (|103p reads: 



tLL 



(n + 1)!' 



9 _ _ 



TV ^ ' ' ' ' (n+ 1)!' 



The evaluation of the second, two-loop graph with the two propagators (|102p in Fig[T] is more 
complicated due to additional divergency of the box subgraph. But technically it can be done in the 
same way therefore we shall not repeat this discussion. 
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